Haldane and dimer gap in general double spin-chain models 
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We have obtained an analytic expression for the k dependence of excitation energy gap for an 
arbitrary double 5* = 1/2 spin-chain by using the nonlocal unitary transformation and the variational 
method. It is checked to explain the gap behavior of various systems, which include the Haldane 
system and the dimer system in both extreme limits, and also the ladder model and the Majumdar- 
Ghosh model. The string order parameter, the dimer order parameter, and the local spin value are 
also calculated in the ground state. The ground-state energy exhibits a great stabilization by an 
antiferromagnetic bond dimerization, which might be realized in various new compounds. We also 
mention the relation of the convergence to the Haldane state with the spin-exchange symmetry of 
the model. The excited state has one domain wall of a local triplet type except in the vicinity of 
the Majumdar-Ghosh point, where a local triplet is decomposed to two S = 1/2 free spins moving 
among the singlet dimers. 
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I. INTRODUCTION 

The low-dimensional quantum systems with the ex- 
citation energy gap have been attracting much interest 
both theoretically and experimentally,^ though the in- 
terest was only from the theoretical side until recently. 
The simplest theoretical spin model may be the dimer 
model that consists of independent pairs of 5* = 1/2 
spins connected by an antiferromagnetic (AF) interaction 
bond. The ground state is a product of a singlet dimer 
state on each bond, and the excitation gap is the singljatp. 
triplet dimer gap. The MaJumdar-Ghosh (M-G) modeflEl 
and the A chain modeloB also realize the perfect singlet 
dimer ground state, and thus the gap is intrinsically the 
dimer gap. Another well-known model that has a differ- 
ent origin of the gapas the S = 1 AF spki|d|iain, so-called 
[ Xhe ladder modelLTEil and the bond 
3~t3 interpolate the dimer model and 
the Haldane system by changing a strength of interac- 
tion bonds as a parameter from +oo to — oo. Therefore, 
these models were investigated mainly to clarify the Hal- 
dane system. Our understandings up to now are that the 
dimer state continuously changes to the S ™ i-Haldane 
state without any explicit phase transitionJrllj On the 
other hand, the phase transition becomes of the first Cf^ 
der in a ladder model with both diagonal interactions .E£l 

Situation has changed since it became possible to syn- 
thesize various cornpouads that actually realize the above 
theoretical models.EjTa For example, magnetic suscepti- 
bility measurements on KCuClsEZl and on CaV205t3 in- 
dicate a spin gap behavior, and the experimental data are 
considered to be explained through the frustrated dou- 
ble spin-chain model. E2I In such an analysis, we need to 
estimate the gap as a function of the strength of interac- 
tion bonds. Then tbe-iSusceptibility can be calculated by 
using the gap valuec^ 



In this paper, we consider the generalized double spin- 
chain system defined by its next-nearest-neighbor inter- 
action, Ji, and the alternating nearest-neighbor interac- 
tion, J2 and J3, as 



N 

Ti- = 2_^ Jli^n • Cri+l + Tri ' '''ri+l) 
n=l 



(1) 



Here, N is the linear size of the system, and \cr\ = \t\ = 
1/2. Figure |] shows the depicted lattice. 
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FIG. 1. Shape of the general double spin-chain model we 
treat in this paper. 

The system is reduced to the M-G model with a choice 
of parameter set {Ji,J2iJ'i) = (0-5,1,1), the isotropic 
ladder model with {Ji,J2,J3) = (1,1,0), or (1,0,1), and 
the S = 1 AF chain in the limit J2 — > —00. The case 
with Ji — corresponds to the bond-alternation modei 
in a single chain. Recently, the string order parametercHI 
and the energy for the ground state of this systeari have 
been calculated by the matrix-product method.c^l Heije, 
we make use of the nonlocal unitary transformationcj E£l 
and give explicit expressions of the ground-state energy, 
the excitation gap, string order parameter, the dimer or- 
der parameter, local spin value in the ground state, and 
the domain wall spin value in the excited state for arbi- 
trary (Ji, J2, Js). This transformation ris an adaptation 
of the Kennedy- Tasaki transformationEJ of the 5 = 1 



system to the double S = 1/2 spin-chain systems, and 
is known to be powerful when the ground state is either 
in the Haldane state or in the state with strong dimer 
correlation. 

In Sec. P, we introduce the transformation and the 
variational method employed in this paper. Then the 
energy, a local bond-spin value and the order parameters 
for the ground state are estimated and compared with 
the numerical diagonalization results of the A^ = 12 lat- 
tice. Section Q] describes the excited states, where we 
consider two types of the variation. One is the local- 
triplet domain wall excitation, which we give an explicit 
form for in the whole phase space. The other one is what 
we call the kink-antikink excitation which is governed by 
two S ~ 1 /? free spins moving among the perfect singlet 
dimers.tStZl This type is the elementary excitation near 
the M-G point. 



II. GROUND STATE 

We first rewrite the Hamiltoniari_Eq. (mj, with the 
nonlocal unitary transformation.EiTO The transforma- 
tion is defined by U in the following. 
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U=Y[Un 



Un = P+ + P- expiiwS^], 
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where P+ (Pn) is the projection operator onto states 
with the even (odd) number of S^ = ±1 for i < n — 1. 
Then the Hamiltonian p) is transformed as. 
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U '^HU = V Ji(-«+i - «+i - 4«+i«+i) 
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^J2^n ■ Tn 



(6) 



We consider the following variational basis for the ground 
state of this Hamiltonian. 
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n— 1 n— 1 

|5n>-(IT,i>-U,T»/V2 (8) 

|r„) - a\ T,T> +/3(l T,i> + I i,T))/V2 + 7l i,i) (9) 

I tiT)'s are the states of |cr^,T^). &, a,/3,7 are the real 
variational parameters and satisfy the normalization con- 
dition, a^ +/3^ +7^ = 1. These parameters are supposed 



to be invariant of n, since we consider the uniform ground 
state. In this sense, the present analysis is variational. 

A state with & = is a singlet dimer state on the <t„- 
Tn bond, a state with h = •\/3/2 is the other singlet dimer 
state on the cr„+i-T„ bond, and a state with & = 1 cor- 
responds to the pure VBS state on the cr„-T„ bond. It 
should be noted that our approximation is not the single- 
site approximation regarding the original Hamiltonian so 
that the (t„_|_i-t„ dimer can be represented by Eq. (|^) 
with h = V3/2. 

The energy expectation value is calculated as 
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+ (2Ji-J3)62(l-62) 

- 3J3&V1 - fc2/3(a^ - 7^) 
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(10) 



We can easily find this minimum value under the con- 
straint, a^ -I- /3^ + 7^ = 1, by using the Lagrange multi- 
plier. The energy expectation value eo is 

eo ^ (j2 - ^foVi) [b' - - ^ J3(6 + ./W^)f 

(11) 
with four possible choices of the parameters (a, (3, 7) as 



(a,/?,7) = (±V273,yT73,0), 
(0,-v/T73,±V273), 

and b determined implicitly through 



J2 



or 




J3 



6 = 0. 
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2Ji (13) 
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The four-fold degeneracy in the choice of (ckjAif 
responds to the degeneracy of the edge states.Ej't 
state with & = is a trivial singlet dimer ground state at 
J2 = 00. The other one, Eq. (13), represents a nontrivial 
state that will be the ground state for the most of the 
parameter space. We solved Eq. (|l^) numerically by the 
bisection method for arbitrary (Ji, J2, J3). 

Before going through the details of the following vari- 
ational results, let us notice that the system possesses 
the symmetry which exchanges J2 bonds and J3 bonds. 
It does not matter if we solve the eigenvalue problem 
exactly, however, the variational results are dependent 
upon this exchange. Therefore, we must do the variation 
on a system whose J2 and J3 values are exchanged as 
well as on a system with a given parameter (Ji, J2, J^), 
and must compare both results. In the present analysis. 



the ground-state energy is always lower if we exchange 
J2 and J3 in the case of J2 > Js- 
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FIG. 2. The ground-state energy for (a) Ji — 0, and J3 — 1 
(bond-alternation model), (b) Ji — 0.5, and J3 = 1 (includes 
the M-G model at J2 — 1), and (c) Ji = 1, and J3 = 1, (in- 
cludes the isotropic ladder model at J2 = 0), respectively. Cir- 
cles denotes the numerical diagonalization results of Af = 12, 
and lines are the variational estimates. 



We consider three cases which only differ a choice of 
Ji. J3 is always set equal to 1, and J2 is a variable that 
we move from -|-cx) to —00. The case (a):Ji = and 
J3 = 1 is the bond-alternation model. It includes a pure 



dimer model at J2 = 0, and the uniform S — 1/2 AF 
spin chain at J2 = 1- The case (b):Ji = 0.5 and J3 = 1 
includes the M-G model at J2 = 1. The case (c):Ji = 1 
and J3 = 1 includes the isotropic ladder model at J2 = 0. 
In addition, these three cases have a pure dimcr model 
at J2 — -foo and the S ^ 1 AF chain at J2 = —00. 
Let us call the singlet dimer state on the J2 bonds the 
J2-singlet, and also that on the J3 bonds the J3-singlet 
hereafter for simplicity. 

Figure ^ shows the J2 dependence of the ground-state 
energy, cq, compared with the numerical diagonalization 
results of a system with N = 12 under the periodic 
boundary conditions. The variational estimates are done 
by exchanging J2 and J3 in the region J2 > J3. The 
energy agrees with the numerical ones fairly well, partic- 
ularly in the dimer region (J2 > J3 = 1). The difference 
becomes visible in the Haldane region (J2 < J3 — 1) as 
Ji increase. The energy takes maximum at the fully- 
frustrated point J2 = J3 = 1 in Fig. g (b) and (c). As 
J2 goes away from J3 = 1, the energy decreases because 
the frustration is relaxed. It should be noted that this 
energy stabilization is stronger in the dimer region com- 
pared with the Haldane region. Thus, we point out here 
a possibility that the ground state of the real double spin- 
chain compounds are also easily stabilized to the dimer 
ground state by a lattice dimerization similar to the spin- 
Peierls system corresponding to the case Fig. 0(a). 

Next, we calculate the local bond-spin value defined by 

(5( J2)) = (cr„ . x„ + ^) = &2 



(5(J3)) = 



fn+l 



-) = --{b + ^3{T^¥~yf 
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Here, S{Ji) denotes the local spin expectation value along 
the Jj-bond. Since a local bond-spin value is not a good 
quantum number, we have to define it by using a pro- 
jection operator which selects local triplet component 
from the composite spin magnitude. The bond corre- 
lation cr„ • T„ + 3/4 serves as this projection operator for 
S{J2)- S{Ji) takes zero for the J^-singlet state, while the 
other bond-spin takes a value of 3/4. 

Figure || shows J2 dependence of S{J2) and S^J^). 
Lines are the variational estimates stated above, and 
symbols are the numerical diagonalization results of iV == 
12. 

Consistency between the variational estimates and 
the numerical results is generally excellent except for 
the SiJs) of Fig. |(c). In Fig. |(a), namely the 
bond-alternation model, there are two trivial pure dimer 
points. The ground state is the J2-singlet state at 
J2 — -l-cxD, and is the J3-singlet state at J2 = 0. These 
two points are equivalent to each other if we exchange the 
J2 bonds and the J3 bonds. Thus, the local bond-spin 
values are symmetric from the isotropic point at J2 = 1, 
where the model reduces to the uniform S' = 1/2 AF spin 
chain. The J3-singlet state continuously changes to the 
Haldane state in the limit of J2 —^ —00, as is visible by 
the S'( J2) converging to 1. In the model that includes the 



(15) 
(16) 



M-G model, Fig. || (b), the ground state is exactly the 
J2-singlet for J2 > 1, where S{J2) = and ^(Ja) = 3/4. 
They show a sudden jump at J2 = 1, since the Js-singlet 
is degenerate with the J2-singlet at this point. Situation 
is rather different in Fig. ra (c), since the Js-singlet never 
becomes the ground state in this case. ^'(Js) increases 
from 3/4 as J2 decrease from +00 until it suddenly de- 
creases at the symmetric point, J2 = J3 = 1- S{J2) now 
takes nearly the triplet value at J2 '--^ 1. At the isotropic 
ladder point of J2 = 0, the diagonal spins of the ladder 
form an almost triplet state since S'( J2) ~ 1, and the rung 
spins are far from the singlet state, since S{J3) ~ 0.3. 



FIG. 3. The local bond-spin value on the J2 bond (S{J2)) 
and that on the J3 bond (^(Ja)) plotted against J2/(l — J2) 
for J2 < 0, and against J2/(l -I- J2) for J2 > 0, where (a) 
Ji = 0, and J3 = 1 (bond-alternation model), (b) Ji — 0.5, 
and Js — 1 (includes the M-G model at J2 = 1), and (c) 
Ji — 1, and J3 = 1, (includes the isotropic ladder model at 
J2 =0), respectively. Lines are the variational estimates. Tri- 
angles and circles are the numerical diagonalization results for 
S'( J2) and S( J3) of a system with 24 spins under the periodic 
boundary conditions. 



We can also estimate tl 
Nijs and Rommelse, Ostr,L 
eter, Odim defined by Hida. 
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In addition, Odim{J2) — 1 only when 6 = and otherwise 
it vanishes; Ostr(^3) = Odim(^2)/4. 

Figure ^ shows the J2 dependence of the dimer and the 
string order parameter in the ground state of three cases 
mentioned above. Symbols denote the numerical diago- 
nalization results of A^ = 12, and lines are the variational 
estimates. The order parameters both on the J2 bonds 
and on the J3 bonds are plotted in the same figure. For 
example, the Ostr(^2) is defined by Eq. (18) when we 
consider the bond-spin Sn is along the J2-bond, and the 
OstriJs) is defined with S'„ along the Js-bond. Hence, 
Ostr{Ji) denotes the den Nijs- Rommelse string order re- 
garding the triplet state of each J^-bond. On the other 
hand, Odim{Ji) expresses the dimer order on each Ji- 
bond. We exchanged J2 and J3 for J2 > J3 as mentioned 
before. 

Consistency of the variational estimates with the nu- 
merical results is excellent except in the limit of J2 -^ 
—00. The variation gives the pure VBS state while the 
numerical ones converge to the correct S — 1 value. This 
is because we used a single-site approximation of the 
transformed Hamiltonian in the variation, and therefore 
our estimates always become worse when the correlation 
length of the ground state is rather long as is the case 
in the Haldane state. Positive values of Odim(J'2) in the 
region J2 < J3 = ^ are the finite-size effect, and should 
vanish in the thermodynamic limit. 

In the system of Fig. H (a), i.e., the bond-alternation 
model, Odim = 1 at the two pure dimer points. The vari- 



ational estimates are particularly good near these two 
points because of the short correlation length, while they 
become poor near the uniform S = 1/2 chain point at 
J2 — 1 and the uniform S = I chain point at J2 = — c». 

In Fig. § (b), the Odim(>/2) is 1 for J2 > 1, since the 
ground state is exactly the J2-singlet. The Ostr(J2) and 
the OdimiJs) is zero in this region. However, the OsniJ^) 
takes a finite value of 1/4. This is a natural consequence 
from the definition of the string order parameter. Thus, it 
is not adequate to discriminate the Haldane phase from 
the dimer phase by vanishing or non-vanishing of this 
parameter alone. We should determine the phase by its 
value ranging from 1/4 in the dimer state to 0.37 in the 
Haldane state. Therefore, it is very difficult to draw a 
phase boundary line in most cases. As J2 decrease from 1, 
the Ja-singlet remains the ground state by changing con- 
tinuously to the Haldane state in the limit of J2 ~* ~oo 
as is observed in the behavior of the OdimiJs) and the 

Ostr(J2). 

The inconsistency between the variational estimates 
and the numerical results in Fig. ^ (b) is larger com- 
pared with that of (a), and becomes distinct in (c). The 
difference in Fig. |4| (c) is already clear in the vicinity of 
J2 = 1, and remains until J2 — > —00. The variation can 
explain the behavior of the order parameters only qual- 
itatively in this plot. For example, Ostr(t^2) is almost 
invariant of J2 in the region J2 < 1. 
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FIG. 4. The Ostr and Odim on the J2 bonds (circles) 
and those on the J3 bonds (triangles) are plotted against 
^2/(1 — J2) for J2 < 0, and against J2/(l + J2) for J2 > 0, 
where (a) Ji = 0, and J3 = 1 (bond-alternation model), (b) 
Ji = 0.5, and J3 = 1 (includes the M-G model at J2 = 1), 
and (c) J\ = 1, and J3 = 1, (includes the isotropic ladder 
model at J2 = 0), respectively. Solid (broken) lines are the 
variational estimate for Ostr (Odim)- 

Looking at Figs. Isl and 0, we notice that the conver- 
gence to the Haldane state becomes faster as Ji increases 
from to 1. For example, the Osti{J2) of Fig. Q (c) takes 
a value of the Haldane state even in the isotropic ladder 
model (J2 = 0). This evidence may allow us to consider 
that the ground state of the isotropic ladder model is 
more like the Haldane state rather than the dimer state. 

We relate this tendency to the symmetry of the model 
with respect to the exchange of two spins that couple to 
form the S = 1 state in the Haldane limit. In general, 
this spin-exchange symmetry is necessary to realize the 
Haldane state, since each 5 = 1 unit should have this 
symmetry. Figure ^ shows a consequence of an exchange 
of the spins <t„ and t„ for three models we have consid- 
ered in this paper. Bold lines denote the J2 bonds, thin 
lines are the bonds of their magnitude 1 , and broken lines 
are those of 0.5. These models only have this symmetry 
in the limit of J2 -^ ±00, and thus the Haldane state 
becomes the exact state only in this limit. However, the 
number of the bonds that are invariant before and after 
this operation increases as Ji increases from to 1. In 
this sense, the case (c) Ji — 1 and J3 = 1 is closer to the 
spin-exchange symmetry. 

On the other hand, the ladder model with both diag- 
onal interactions as depicted in Fig. @(d) has this sym- 
metry for arbitrary J2 values. As shown in Fig. |^(e), the 
first-order transition from the dimer. phase to the Hal- 
dane phase occurs at J2 — 1.40148.E3 In this figure, the 
symbols denote the numerical results of a system with 20 
spins. Negative values of the dimer order are the finite- 
size effect. Therefore, we conclude that the convergence 
to the Haldane state becomes faster as the system gains 
this spin-exchange symmetry. At the same time, our vari- 
ational estimates fails faster. 
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FIG. 5. Symmetry of the model with respect to the 
spin-exchange of (t„ and t„ for (a) Ji — 0, J3 = 1, 
(b)Ji = 0.5, J3 = 1, (c)Ji = 1,J3 = 1, and (d) the lad- 
der model with both diagonal interaction bonds. Bold lines 
denote the J2-bond, thin lines are the bonds with their mag- 
nitude 1, and broken lines are those of 0.5. (e)The string and 
the dimer order parameter of the model (d). 



III. EXCITED STATE 

The elementary excitation in one dimension is intrinsi- 
cally a state with one domain wall between the degener- 
ate ground states. Within the present variational scheme, 
we consider the following one domain- wall state under the 
open boundary conditions. 

/ i N \ 



i \n— 1 n— i+1 / 

= J2C^^^■ (19) 



Here, the parameter sets (a,/3, 7) and (a',/?',7')are any 
two of the four possible choices given in Eq. (|l^), and 
b is determined by Eq. (p^). For example, we use 
the set (a,/3,7) = {y/2/3, ^/T/3,0) and (a',/3',7') = 
(— 1/2/3, 1/1/3,0). Of course, the choice does not af- 
fect the final results. A domain wall is located between 
the ith site and the (i + l)th site. This definition of the 
trial function becomes equivalent to the solitonic jOxcita- 
tion of Fath and SolyornEl in the AKLT model.y The 
spin expectation of the domain wall is defined by 



{■)pi\Ti ■ (Ti+l 
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3/4^',) 
(6 + V3(l-62))' ^ 



rOd 



(20) 



Figure g shows the J2 dependence of this value for three 
different cases discussed in the previous section. 
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FIG. 6. The variational estimates of the spin expectation 
of a domain wall in the excited state are plotted against 
J2/(l — J2) for J2 < 0, and against J2/(l + J2) for J2 > 0, 
where (a) Ji = 0, and J3 = 1 (bond-alternation model), (b) 
Ji = 0.5, and J3 = 1 (includes the M-G model at J2 — 1), 
and (c) Ji — 1, and J3 = 1, (includes the isotropic ladder 
model at J2 =0), respectively. 

The excitation becomes a local triplet at the domain 
wall when the ground state is exactly the singlet dimer 
state. As the ground state changes to the Haldane state, 
the local triplet smears out and consequently the local 
spin value at the domain wall decreases, since the total 
spin of the excited state is always 1 in this system. In 
the Haldane limit, it takes a value of 31/36 ^ 0.861. 

The basis relations and the matrix element of the 
Hamiltonian is calculated as 
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The variation (*i|7Y|*i)/(*i|*i) can be calculated 
by the Fourier transformation, since the denominator 
(^il^i) is diagonalized in the thermodynamic limit. 
Namely, by \4>k) = Z^n *3xp[ifcn]|^„). Then the energy 
gap is obtained with respect to the wave number of the 
domain wall k as, 
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where a is defined by Eq. (pi|). 

We plot this estimate with k ~ and fc = tt in Fig. 
and compare with the numerical results. 
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7. J2 dependence of the energy gap obtained by the 
numerical diagonalization (symbols) , and the variation (lines) 
plotted against J2/(l — J2) for J2 < 0, and against J2/(l + J2) 
for J2 > 0, where (a) Ji — 0, and J3 = 1 (bond-alternation 
model), (b) Ji = 0.5, and J3 = 1 (includes the M-G model at 
J2 = 1), and (c) Ji = 1, and J3 — 1, (includes the isotropic 
ladder model at J2 =0), respectively. We also plot the esti- 
mate by the kink-antikink variation (broken dash line) in (b) . 



In the vicinity of the M-G point, {Ji,J2,J3) — 
(0.5, 1, 1), the lowest excitation is a kink-antikink state, 
which consists of {N — 1) singlet dimer pairs and two free 
5=1/2 spins .B We call these free spins a kink and an 
antikink. In such a case, we try another variation. These 
two free spins are mobile in general and thus we have to 
consider the following trial function. 



j-i 



N 



iV'y) = n I'^w) n !"("' ^^ ^' ^)) n i^w)' (2?) 



under the periodic boundary conditions. Here, |n.(0)) 
stands for the singlet dimer state on the CTn-Tn bond 
with 6 = 0. |rt(a, /3, 7, b)) becomes another singlet dimer 
state on the cr„+i--r„ bond when h = V3/2 at the M-G 
point. Therefore a kink or an antikink should exist at 
the domain wall. 

In the A chain, a kink becomes localized and thus 
the problem can be reduced to a one-body problem of 
a moving antikink. We can solve it exactly-|paiid the 
wave function is given by the Airy function.r^rl Here, 
we solved this variation only numerically for a finite sys- 
tem of TV = 30. The results are also plotted in Fig. CTb) 
by a broken dash line. 

The relative motion of this two-body problem is same 
as the one-body problem of the A chain. That is, we 
obtain the same equation if we set one free spin at the 
origin and rewrite the variational problem with respect 
to the one body-problem of the other free spin. This is 
also verified by that the wave function obtained numeri- 
cally is well fitted by the Airy function, and that the gap 
enhancement, £'gap( J2) — £^gap(^2 = 1), is almost equal to 
that of the A chain; it obeys a power law with exponent 



2/3.c3EZll23 Of course, the gap itself cannot be obtained 
correctly only by the relative motion but together with 
the motion of the center of mass. 

Figure shows the J2 dependence of the energy gap 
estimated above for the cases with (a) Ji = 0, J3 = 1, (b) 
Ji — 0.5, J3 = 1 and (c) Ji = 1, J3 = 1. The local triplet 
excitation with k — t: is depicted by solid lines, that with 
fc = is by broken lines, and the kink-antikink excitation 
is by a broken dash line in Fig. R(b). The lowest gap in 
the k = TT sector and that in the fc = sector calcu- 
lated by the numerical diagonalization of an A^ = 12 lat- 
tice are depicted by circles and by triangles, respectively. 
The variational estimates with k = n are quite excel- 
lent for all the plots. We consider this is because the 
local approximation is usually good for the wave num- 
ber TT, which changes the phase of the wave function by 
only one lattice spacing. Only a difference is that the 
local triplet variation converges to the VBS value in the 
J2 — *■ —00 limit, while the numerical one converges to the 
Haldane value. When the ground state is exact singlet 
dimer state (J2 > J3 = 1 in Fig. 0(b)), our variational 
excitation is a local triplet state and thus the gap value 
is equal to J2 with no dispersion. On the other hand, 
the estimates with A; = only explain the gap behavior 
qualitatively, and become worse with an increase of Ji. 
They are only valid near the pure dimer points in Fig. 

a) . In the vicinity of the gapless point of J2 = J3 = 1 
(S=l/2 chain), our variational scheme breaks down, and 
thus the gap estimate should be done by aiiOLLher method. 
For example, the bosonization techniqueLj Ej gives 



0( 



18 



|2x-l| 
l-x 



(28) 



with X = >/2/(l + J2), and agrees with our numerical re- 
sults very well. This is plotted by a broken dash line in 
Fig. 0(a). The variational estimate of the kink-antikink 
type is also consistent with the numerical results. We 
can consider that it is valid as long as the singlet dimer 
state is exactly or approximately the ground state. 



amount of this symmetry, the convergence of the ground 
state to the Haldane state becomes faster, or in other 
words, the transition becomes close to the first-order. In 
the case of the system with^full symmetry, the transition 
is strictly of the first-order .113 

The excited state is formulated by a domain wall be- 
tween two of the four-fold degenerate ground states. In 
the dimer region, this domain wall is equivalent to a local 
triplet. We obtained an explicit form for the dispersion 
relation of the gap for arbitrary Ji, J2, and J3. That 
is, we first calculate the value of b by Eq. ( p^ ) for a 
given (Ji, J2, J3). Then, the dispersion is given by Eq. 
(EQ) with a = 4/3&^ — 1. We confirmed that our vari- 
ational estimate is generally good in the dimer region, 
J2 > J3, and is especially excellent for the excited state 
with k = TT. The lowest excitation ncar|t|ba|M-G point is 
well-explained by a kink-antikink state.Hc3'E] 

Our variation employs the single-site approximation in 
the transformed system, which of course becomes worse 
when the ground state has rather long correlation length, 
e.g. like in the Haldane state. Therefore, our estimate 
fails faster with the convergence to the Haldane system. 
We must go beyond the single-site approximation for the 
sake of quantitative agreements in this region. 

Finally, we point out a possibility that the ground state 
of a real compound is stabilized by a lattice dimerization 
to the dimer state, since its energy stabilization is quite 
significant as we observed in Fig. g. In fact, the suscep- 
tibility of KCuCls is quantitatiiiely explained by a single 
dimer model with J =:r43.8K.E3 The situation is similar 
in the case of CaV20E 
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IV. SUMMARY AND DISCUSSION 



We have investigated the general S = 1/2 double spin 
chain (J1-J2-J3 model) by means of the nonlocal unitary 
transformation and the variation. The model includes 
the dimer model and the Haldane system in its both ex- 
tremes. A ground-state change occurs at the symmetric 
point, J2 = J3 = 1, as is visible in the string and the 
dimer order parameter. The ground state of J2 < J3 con- 
tinuously changes to the Haldane state in the J2 — > — 00 
limit, and the other ground state of J2 > J3 becomes 
the dimer state in the J2 ^ 00 limit. Two states are 
degenerate at the symmetric point. 

We relate the convergence to the Haldane state with 
the symmetry of the exchange of spins that will couple to 
the S = 1 state in the Haldane limit. As the system gains 
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